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Abstract. Computing a polygon defining a set of planar points is a classical problem of 
modern computational geometry. In laboratory experiments we demonstrate that a concave 
hull, a connected a- shape without holes, of a finite planar set is approximated by slime 
mould Physarum polycephalum. We represent planar points with sources of long-distance 
attractants and short-distance repellents and inoculate a piece of Plasmodium outside the 
data set. The Plasmodium moves towards the data and envelops it by pronounced proto- 
plasmic tubes. 
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1. Introduction 

A slime mould Physarum polycephalum has a rich life cycle |30|, which includes fruit 
bodies, spores, single-cell amoebas, and syncytium. Plasmodium is a vegetative stage 
of Physarum polycephalum, syncytium, a single cell with many nuclei. The Plasmodium 
consumes microscopic particles. During its foraging behaviour the Plasmodium spans scat- 
tered sources of nutrients with a network of protoplasmic tubes. The protoplasmic network 
is optimised to cover all sources of food and to deliver robust and speedy transportation 
of nutrients and metabolites in the Plasmodium body. Plasmodium's foraging behaviour 
can be interpreted as computation, when data are represented by spatial configurations of 
attractants and repellents, and results by structures of protoplasmic network |3, 2|. Plas- 
modium satisfactory solves many computational problems with natural parallelism, includ- 
ing shortest path [^^^l. implementation of storage modification machines |4 |, Voronoi 
diagram logical computing (331 0, process algebra |28|; see overview in |2|. 

A convex hull of a finite planar set is a smallest region of the plane containing the 
set l27ll . Classical algorithms of convex hull construction include Jarvis' gift wrapping flTl , 
Preparata-Hong's divide and conquer |26|, Graham's scan of pre-sorted set of points flgl, 
and Akl-Toussaint's quick algorithm based on removing points lying inside a convex quadri- 
lateral of extreme points Q ; see reviews and discussion of the algorithms in (27] [TTl 
and codes of implementations in |25|. Multi-processor and systolic-processor solutions 
of the convex hull problem are mostly based on parallelisation of the classical serial al- 
gorithms ilOl El or exploring extremal properties of the given set flSl EH [121 ; see also 
review in fS^. Cellular- automata algorithm are based on propagating patterns halting their 
growth in the hull's boundaries 1 1 1 or encapsulating the data set and developing synchroni- 
sation patterns |^30|, are theoretical precursors for bio-inspired approach to computation 
of convex hulls. 

A convex hull is economical but rough reconstruction of an object shape form the ob- 
ject's sample points. It often lacks details necessary for realistic visualisation of objects. 
In 1983 Edelsbrunner, Kirkpatrick and Seidel 1 13J introduced a-shapes, generalisation of 
convex hulls, where control parameter a allows for a smooth transition between crude 
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and fine approximations of point sets. There are hundreds of algorithms for computing 
a-shapes yet no experimental laboratory prototypes of chemical, physical or biological 
computers for approximation of a-shapes. Present paper partly fills the gap and shows 
how to approximate a-shapes without holes using foraging behaviour of R polycephalum. 

2. Experimental 

A convex hull of a finite set P of planar points is the smallest convex polygon that 
contains all points of P (Fig. [T^). a-hull of P is an intersection of the complement of 
all closed discs of radius 1/a that includes no points of P a-shape is a convex 

hull when a ^ oo. With decrease of a the shapes may shrink, develop holes and become 
disconnected, the shapes collapse to P when a ^ 0. A concave hull is non-convex polygon 
representing area occupied by P. A concave hull is a connected a-shape without holes 
(Fig.[T]3). 

Problem 1. Given planar set P represented by physical objects Plasmodium of P. poly- 
cephalum must represent concave hull ofV by its largest protoplasmic tube. 

We cultivate P. polycephalum in plastic containers, on paper kitchen towels sprinkled 
with still drinking water and fed with oat flakes. For experiments we use polystyrene Petri 
dishes (round, diameter 120 mm, and rectangular 120 x 120 mm) and 2% agar gel (Select 
agar, Sigma Aldrich) as a non-nutrient substrate. Images of Plasmodium are recorded 
by scanning Petri dishes in Epson Perfection 4490. Photos are taken using FujiPix 6000 
camera. 

3. Results 

The fist algorithm of convex hull construction 1 17 | was based on cognitive tactic tech- 
niques we use in our everyday 's life. We select a starting point which is extremal point of 
P. We pull a rope (anti-)clockwise to other extremal point. We continue until the set P is 
wrapped completely. The computation stops when we reach the starting point. Let we rep- 
resent data points P by sources of attractants only, e.g. by oat flakes (Fig.[T]:de). We place 
a piece of Plasmodium at some distance away from the set of points, see scheme in Fig.[T^. 
The Plasmodium propagates towards set P, colonises oat flakes (Fig. [T]:d) and spans them 
with a network of protoplasmic tubes (Fig. [T^). No hull is constructed. When all data 
points P are colonised and spanned by protoplasmic network the Plasmodium ventures to 
explore the space around P. 

Proposition 1. P. polycephalum does not compute concave or convex hull of a set repre- 
sented by attracting sources. 

Using repellents only, e.g. as a set of obstacles between attractants and inoculation 
site, will do no good because Plasmodium will just pass around the repellents and leave 
(see details in |2|, controlling Physarum with salt). The only solution would be to em- 
ploy attractants to 'pull' Plasmodium towards planar set P and to use repellents to prevent 
Plasmodium from spanning the points of P (Fig. [2^). Strength of repellents should be 
proportional to a and thus will determine exact shape of the constructed hull. This corre- 
sponds to original definition 1 13 1 that a-hull of P is the intersection of all closed discs with 
radius / / a that contain all the points of P. 

Proposition 2. Plasmodium of P. polycephalum approximates connected a-hull (without 
holes) of a finite planar set, which points are represented by sources of long-distance at- 
tractants and short-distance repellents. 
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Figure 1 . Basics, (a) Example of convex hull, (b) Example of concave 
hull. Points of set P are empty discs, (c-d) Spanning oat flakes, repre- 
senting P, by a network of protoplasmic tubes: scanned image (c) and its 
binarisation (d). (e) Scheme of the Plasmodium propagation, circles are 
points, or oat flakes, of P and star marks site of inoculation, arrows are 
protoplasmic tubes. 
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Figure 2. Experimental, (a) Proposed distribution of attracting and 
repelling gradients which may force Plasmodium to approximate a con- 
cave hull (arrows aiming towards set of discs are attractive forces, arrows 
originating in data points, discs, are repelling forces), (be) Plasmod- 
ium's interaction with the single half-pill of Kalms Tablets: scanned (b) 
and binarised (c) images; at the beginning of experiment Plasmodium 
is inoculated in southmost part of Petri dish (pill of Kalms Tablets is at 
the centre of the Petri dish), (d) Typical experimental setup, pills repre- 
senting P and the initial site of Plasmodium's inoculation are shown by 
arrows. 
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To prove the proposition constructively we must find an appropriate representation of 
P and demonstrate in experiments viabiHty of the approach. While experimenting with 
different potential candidates for combined representation of attractants and repellents we 
found that the Plasmodium's reaction to pills of Kalms Tablets and Kalms Sleep[^is some- 
what unusual. When presented with a half-pill of the Kalms Tablets/Sleep the Plasmodium 
propagates towards the pill and forms, with its protoplasmic tubes, a circular enclosure 
around the pill (Fig.[2|)c). Such a unique behaviour of Plasmodium in presence of Kalms 
Tablets/Sleep indicates that a Plasmodium could implement Jarvis's Gift Wrapping algo- 
rithm ITtI, adapted to concave hulls without holes, if points of P are represented by the 
pills. 

We tested feasibility of the idea in 25 experiments. All experiments were successful. 
In each experiment we arranged 4-8 half-pills (they represented given planar set P) in a 
random fashion near centre of a Petri dish and inoculated an oat flake colonised by Plas- 
modium 2-4 cm away from the set P. A typical experiment is illustrated in Fig. [3^-d. In 
12 h after inoculation plasmodium propagates towards set P and starts enveloping the set 
with its body and network of protoplasmic tubes (Fig.[3^b). The plasmodium completes 
approximation of a shape by entirely enveloping P in next 12 h (Fig.[3]:d). The plasmodium 
does not propagate inside configuration of pills (Fig. [3^). 

Configuration of P in Fig. |3^-d favours approximation of a convex hull. If spatial 
configuration of points curves inwards then concave hull is approximated (Fig.|4]). Typical 
hulls approximated by plasmodium of P. polycephalum in our experiments are shown in 
Fig.|5] Hulls constructed by plasmodium match their counterparts calculated by classical 
algorithms. 

Let us verify experimental laboratory results with computer simulation. A profile of 
Plasmodium's propagating zone is isomorphic to shapes of wave-fragments in sub-excitable 
media |2 |. When active zone of P. polycephalum propagates two processes occur simulta- 
neously — propagation of the wave-shaped tip of the pseudopodium and formation of the 
trail of protoplasmic tubes. We simulate the (chemo-)tactic traveling of plasmodium using 
two-variable Oregonator equation (161 Ell • 

ot £ u-\-q 

dv 
dt 

The variable u is abstracted as a local density of Plasmodium's protoplasm and v reflects lo- 
cal concentration of metabolites and nutrients. We integrate the system using Euler method 
with five-node Laplasian operator, time step = 5 • 10~^ and grid point spacing Ax = 0.25, 
with the following parameters: = 0o - ^/2, A = 0.0011109, 0o = 0.0766, £ = 0.03, 
/= 1.4,^ = 0.022. 

Parameters q and / are inherited from model of Belousov-Zhabotinsky medium, (j) is 
proportional to local concentration of attractants and repellents. The parameter r] corre- 
sponds to a gradient of chemo- attractants emitted by data planar points. The medium is 
perturbed by an initial excitation, where a 1 1 x 1 1 sites are assigned u=l.O each. The per- 
turbation generates a wave-fragment alike propagation of protoplasm which travels along 
gradient T] . Repellents emitted by sites of P show very limited diffusion. Therefore repel- 
lents can be regarded as impassable obstacles (Fig. [6^). 
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Figure 3. Computation of concave hulls of point set P. Images (ac) 
and their binarizations (bd) taken 12 h (ab) and 24 h (cd) hours after in- 
oculation, (e) Scheme of the Plasmodium interaction with P: points/pills 
of P are circles; diffusing repellents, visible as black halo's surround- 
ing pills in (a) and (c), are shown by grey gradients; feeding boundary 
of protoplasm is marked by dotted line; major protoplasmic tube, which 
approximates C//(P), is shown by solid line. 
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Figure 4. Binarized image of Plasmodium approximating concave hull 
of P. 




Figure 5 . Contours of convex and concave hulls constructed on small 
sets of planar points in laboratory experiments. 

To imitate formation of the protoplasmic tubes we store values of u in matrix L, which 
is processed at the end of simulation. For any site x and time step ^ if > 0. 1 and Lx = 
then Lx = I. The matrix L represents time lapse superposition of propagating wave-fronts 
(Fig. [6|3). The simulation is considered completed when propagating pattern envelops P 
and halts any further motion. At the end of simulation we repeatedly apply the erosion 
operation |2 | (which represents a stretch- activation effect |20| necessary for formation of 
Plasmodium tubes) to L. The resultant protoplasmic network (Fig. [6]:) provides a good 
phenomenological match for networks recorded in laboratory experiments. 
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Figure 6. Simulating Plasmodium approximation of convex hull of P 
in Oregonator model, (a) Gradients of repellents generated by points 
of P. (b) Time lapsed images of propagating plasmodium patterns, 
(c) Structure of protoplasmic tubes developed. 

4. Discussion 

In last five years we witnessed a substantial progress in design and fabrication of work- 
ing prototypes of amorphous biological computing devices based on foraging behaviour 
of slime mould P. polycephalum. Experimental laboratory prototypes of Physarum-based 
processors are developed for computation of minimum spanning tree, relative neighbour- 
hood graph, Delaunay triangulation, and Voronoi diagram, see overview in |f2 |. By present 
experimental demonstration of planar shape (convex and concave hulls) computing by P. 
polycephalum we draw a logic conclusion of the years of studies into computational prop- 
erties of the slime mould with regards to problems of computational geometry. 
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We used Kalms Tablet and Kalms Sleep pills (both brands cause similar effects to Plas- 
modium behaviour) to represent data points. We exploited two properties of the pills. First, 
the pills contain Valerian and hops extracts, which are long distant attractants for P. poly- 
cephalum |6| and thus 'pull' Plasmodium towards the date set. Second, the pills contain 
short-distance, slow diffusing components (potential candidates are magnesium stearate, 
stearid acid and titanium dioxide present, and more likely sucrose [35 1) which may play a 
role of repellents preventing plasmodium from spreading inside the data set. Combination 
of these two properties causes P. polycephalum to develop a network of protoplasmic tubes 
with major tube approximating convex hull. Exact mechanism of short-distance repelling 
is unclear and will be a topic of further investigations. 

Plasmodium of P. polycephalum approximates shape of P by propagating simultane- 
ously clock- and anti-clockwise (two branches). These branches fuse when meet one an- 
other. Thus the slime mould approximates shape of P in time n/2, where n is a perimeter 
of the shape. The time can be reduced further by n by inoculating plasmodium in n loci 
of space at the same. Indeed the slime mould based experimental computation of planar 
shape can not and will not compete with existing algorithms, however will contribute to- 
wards design of future parallel embedded processors made of non-linear chemical media, 
and also in control and navigation of amoeboid robots. 
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